Dynamical growth of the hadron bubbles during the quark-hadron phase transition 
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The rate of dynamical growth of the hadron bubbles in a supercooled baryon free quark-gluon 
plasma, is evaluated by solving the equations of relativistic fluid dynamics in all regions. For a 
non- viscous plasma, this dynamical growth rate is found to depend only on the range of correlation 
£ of order parameter fluctuation, and the radius R of the critical hadron bubble, the two length 
scales relevant for the description of the critical phenomena. Further, it is shown that the dynamical 
prefactor acquires an additive component when the medium becomes viscous. Interestingly, under 
certain reasonable assumption for the velocity of the sound in the medium around the saddle config- 
uration, the viscous and the non-viscous parts of the prefactor are found to be similar to the results 
obtained by Csernai-Kapusta and Ruggeri- Friedman (for the case of zero viscosity) respectively. 
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I. INTRODUCTION 

The phenomena of phase transition has attracted many 
researchers from diverse areas due to many interest- 
ing and common features that occur near the transition 
point. Recently, a considerable amount of attention is 
being paid to the study of relativistic heavy ion colli- 
sions where a phase transition is expected from the nor- 
mal nuclear matter to a deconfined state of quarks and 
gluons jy. The quark gluon plasma (QGP), if formed, 
would expand hydrodynamically and would cool down 
until it reaches a critical temperature T c where a phase 
transition to hadron phase begins. Although the order 
of such a phase transition remains an unsettled issue, a 
considerable amount of work has been carried out to un- 
derstand the dynamics assuming it to be of first order 
and also assuming that the homogeneous nucleation is 
applicable J^-§]. In the ideal Maxwell construction, the 
temperature of the plasma remains fixed at T c during 
the phase transition until the hadronization gets com- 
pleted. However, if the hadronization proceeds through 
nucleation, it will not begin at T — T c due to the large 
nucleation barrier. The nucleation of the hadron bubbles 
can begin only from a supercooled metastable state. If 
the amount of supercooling is small, the nucleation rate 
H is computed from I = A exp (-AF/T) which gives 
the probability per unit time per unit volume to nucleate 
a region of the stable phase (the hadron phase) within 
the metastable phase (the QGP phase). Here AF is the 
minimum energy needed to create a critical bubble and 
the prefactor A is the product of statistical and dynam- 
ical factors. The statistical factor fio is a measure of 
both the available phase space as the system goes over 
the saddle and of the statistical fluctuations at the saddle 
relative to the equilibrium states. The dynamical prefac- 
tor k gives the exponential growth rate of the bubble or 
droplet sitting on the saddle. 

In an earlier work, Langer and Turski || derived the 
dynamical growth rate (ft) of the liquid droplet based on 
a non-relativistic formalism. Subsequently, re was derived 
both by Turski-Langer M and Kawasaki |l for a liquid- 



gas phase transition near the critical point, to be 
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which involves the thermal conductivity A, the surface 
free energy a, the latent heat per molecule I and the 
density of the molecules in the liquid phase n/. The in- 
teresting physics in this expression is the thermal con- 
ductivity which appears as an essential ingredient for the 
transportation of the latent heat away from the surface 
region so that the droplet can grow. For a relativistic 
system like a baryon free quark gluon plasma which has 
no net conserved charge, the thermal conductivity van- 
ishes. Hence, the above formula obviously can not be 
applied to such systems. Therefore, Csernai and Ka- 
pusta re-derived re for a baryon free plasma using earlier 
formalism of Langer- Turski Q , but extending their work 
to the relativistic domain |^| . In the work of Langer and 
Turski ||, re was derived by solving a set of linearized 
hydrodynamic equations in the liquid, vapor and the in- 
terfacial regions. However, in the relativistic formalism, 
Csernai-Kapusta mostly concentrated in the interfacial 
region. Their primary motivation was to know the veloc- 
ity profile in the surface region which was then used to 
estimate the energy flow across the surface. Then they 
used the condition, that the energy flux which is to be 
transported outwards should be balanced by the viscous 
heat dissipation as follows 



dR 

Alu — = -(477/3 + C) vdv/dr, 
dt 



(2) 



where R is the radius of the hadron bubble and v(r) is 
the flow velocity just outside the surface of the bubble. 
Accordingly, they obtained an expression for re given by 
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where rj and £ are the shear and bulk viscosity coefficients 
respectively and Alu is the difference in the enthalpy den- 
sities of the plasma and the hadronic phases, lu = e + p. 
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The above expression implies that energy flow wv is pro- 
vided by the viscous effects. There will be no bubble 
growth in the case of an ideal plasma with zero viscos- 
ity. Thus, the viscosity plays the same role as thermal 
conductivity in case of a relativistic fluid like the quark 
gluon plasma with zero baryon density. This approach 
has been extended by Venugopalan and Vischer jlO| for 
the case of baryon rich QGP where both viscous damping 
and thermal dissipation are significant. On the contrary, 
Ruggeri and Friedman (RF) argued that the energy 
flow does not vanish in absence of any heat conduction or 
viscous damping. Since the change of the energy density 
e in time is given in the low velocity limit by the conser- 
vation equation ||, de/dt — -V.(wv) which implies that 
the energy flow cx wv is always present. Therefore, fol- 
lowing a different approach, Ruggeri and Friedman ]Tl| ] 
derived an expression for k which does not vanish in the 
absence of viscosity, given by 
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The viscous effects cause only small perturbations to the 
above equation. This result is in contradiction with the 
expression given in Eq. (^) according to which the hadron 
bubble will not grow in the limit of vanishing viscos- 
ity. The difference between Csernai-Kapusta (CK) and 
Ruggeri-Friedman (RF) results are due to the technical 
differences in the treatment of the pressure gradients and 
it needs further investigation. Motivated by this, we re- 
derive k using Csernai-Kapusta formalism which is a rel- 
ativistic generalization of Langer-Turski (LT) procedure 
|Q| . However, unlike Csernai-Kapusta, we solve the lin- 
earized hydrodynamic equations in all regions namely the 
exterior quark region, the interior hadron region as well 
as the interfacial or the surface region. We found that in 
the limit of zero viscosity, our prefactor k depends only 
on two scale parameters, the correlation length £ and the 
critical radius of the hadron bubble R. We have also ob- 
tained the prefactor for a viscous medium where it can 
be written in a simple way as the sum of a viscous and a 
non-viscous terms. Interestingly, using certain assump- 
tion for the velocity of sound in the medium around the 
saddle configuration, the viscous and non- viscous compo- 
nents are found to be similar to the results as obtained 



by Csernai-Kapusta 
zero viscosity, Eq. (P 



Eq. (J3|)] and Ruggeri-Friedman [for 
J)] respectively. 
The paper is organized as follows. We begin with a 
brief review of the Csernai-Kapusta and Turski-Langer 
formalism describing the energy-momentum conserving 
equations of motion in section II. In section III, we solve 
these equations to derive the dynamical prefactor. Fi- 
nally, the numerical calculations and the conclusions are 
presented in sections IV and V respectively. 



II. THE RELATIVISTIC HYDRODYNAMICS 
FOR BARYON FREE PLASMA 

In the case of relativistic hydrodynamics, we consider 
the energy density e(r, t) and the flow velocity v(r, t) of 
the fluid as two independent variables that describe the 
dynamics of the system. The equations of motion can be 
obtained from the local conservation laws: 



3^ v = 9 M n" = 0. 



(5) 



Here T^ v is the energy momentum tensor and repre- 
sents the baryon four vector. In the presence of viscosity, 
the energy-momentum tensor T^ v and baryon four vec- 
tor current n M can be decomposed into an ideal and a 
viscous part |l3|JT^ 1 



= [(e + p)u»u v - pg^) + - 



(G) 



(7) 



Here e, p and n are the energy density, pressure and 
particle number density. The fluid four velocity is given 
by u M = 7(1, v) and r Miy and are the dissipative cor- 
rections. The form of the dissipative terms t^ v and 
depend on the definition of what constitutes the local 
rest frame of the fluid. The four velocity should be 
defined in such a way that in a proper frame of any given 
fluid element, the energy and the number densities are 
expressible in terms of other thermodynamic quantities 
by the same formulae, when dissipative processes are not 
present. It is also necessary to specify whether u M is the 
velocity of energy transport or particle transport. Ac- 
cordingly, there exist two definitions for the rest frame; 
one due to Landau and other due to Eckart. In the Lan- 
dau approach, m m is taken as the velocity of energy trans- 
port so that energy three flux T 0i vanishes in a comov- 
ing frame |L^[L4j. In the Eckart definition, it M is taken 
as the velocity of the particle transport and the particle 
three current, rather than the energy three flux vanishes 
in the fluid rest frame So in the Eckart definition 

of rest frame, the particle four vector can be written as 
?i M = (n, 0), whereas in the Landau definition of rest 
frame = (n, v). Therefore, the two frames are related 
by a Lorentz transformation with a boost velocity u/n. 
It is found that due to ill defined boost velocity [[15| , the 
energy three flux in the Eckart frame (which involves heat 
conductivity A) is not well defined as A diverges in the 
limit of chemical potential fi — >0. On the other hand, in 
the Landau definition heat conduction enters as a cor- 
rection to baryon flux. It was shown that inspite of the 
divergence of A, the correction to the baryon flux is 
finite [|l5| . Therefore, we will use the Landau defini- 
tion for the subsequent study and also we will assume a 
baryon free plasma for simplicity. We can now write the 
equations of motion from the conservation law d^T^ = 
using Landau definition pjjji 
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<9 t e = -V.(wv) + 0(v 2 ), 



d t (uv) 



-V.(wv (g) v' 

'4 
-, T 1 + 



Vp 
C) V.v 



0(v 2 



(8) 



(9) 



Here to = (e + p) and rj and £ are the shear and bulk 
viscosity coefficients respectively. We have also assumed 
the low speed limit where 7 f» 1. Although, the fluid 
velocity is small, the velocity of individual particles is 
large. Thus, the expressions for the energy density, pres- 
sure etc. are taken to be same as that in the relativistic 
case. 

Following Ref. @, Eq. © can also be written in terms 
of the Helmholtz free energy /(e) and the usual gradient 
energy iif(Ve) 2 as 

d t (wv) = -V.(wv ® v) - Vp' 

'4 



+ V o ? / + C V.v + 0(v 2 ) 



(10) 



where 



Vp' = -K(V 2 e)Ve 



df 
de 



Ve. 



(11) 



The constant K is related to the surface tension a as 



K 



00 , (de 
dr , , 
ar 



(12) 



It can be noted by comparing Eq. (10) with Eq. (||) 
that — Vp' is not simply a pressure but a combination of 
V/ and a force term — if(V 2 e)Ve which is related to the 
surface tension given by Eq. (|l^). The pressure inside the 
interface differs from that outside so there is necessarily 
a pressure gradient at the interface. The term given by 
Eq. (^Tj) is needed to balance the differential pressure 
otherwise the Euler or the Navier Stokes equation would 
require a changing fluid velocity even in a stationary con- 
figuration, which is unphysical. 



III. SOLUTION OF THE RELATIVISTIC 
HYDRODYNAMIC EQUATIONS 

The above hydrodynamic equations [Eq. (|^) or Eq. 
([nf)] can be solved after linearizing around the saddle 
configuration. The saddle point corresponds to the sta- 
tionary solution when e(r, t) = e(r) and v(r, t) = and 
also e satisfies, M 



K(W 2 e) 



de 



0. 



(13) 



We can now write the equations of motion for small de- 
viations about the stationary configuration by defining 
e = e(r) + v(r, t) where v is a small fluctuation in energy 



density and v = + v(r, f) and linearizing Eqs. (||) and 
( |To| ) around this configuration 



fti/(M) = -V.(wv(r,t)), 
d t (cuv(r, *)) = Ve (-KW 2 + /") i/(r, t) 



+ ( V.v(r,i) 



(14) 



(15) 



Here /" = d 2 f/de 2 , evaluated around the stationary con- 
figuration. The dynamical prefactor k is determined with 
the radial perturbations of the form 



i/(r,t) = v{r)e Kt , 
v(r,t) = v(r)e Kt . 



(16) 



It can be seen from Eqs. ( |14| ) and (15) that the radial 
deviations are governed by the equations of motion of 
the following form 



nv(r) — —V. (wv(r)) , 
(tiiv(r) = Ve (-KV 2 + /") i/(r) 



+V [ ( + V.v(r) 



(17) 



Following Langer and Turski, we find the solution for 
v{r) in each of the three regions (i) the interior region 
of hadron phase, r < R — £ (ii) the exterior region of 
QGP phase, r > R + £ and (iii) the interface region, 
_R — ^<r<i? + ^, where R is the radius of the hadron 
bubble with origin at r = 0. The interfacial region has a 
thickness of the order of the correlation length £. Further, 
it is assumed that, everywhere outside the droplet, the 
energy density e(r) has the value e q , the quark density. 
Within the droplet, e(r) is equal to the hadron density e^. 
Thus, e(r) describes a smooth interfacial profile at r = R 
going from e^ to e q within a region of roughly of the order 
of the correlation length £. We then evaluate the relative 
amplitudes in the above three regions by matching the 
values at the boundaries. Finally, we evaluate n applying 
the condition 



r 2 v{r)dr = 0. 



(18) 



which is the conservation law implied by Eq. (|lj) 

Before we proceed further, it may be noted here that 
the above set of linear equations [Eqs. (|l7|)] is obtained 
from Eq. (10) which contains /(e) and K explicitly. As 
will be shown subsequently, this form is suitable for the 
interfacial region where Ve is nonzero. Such an equation 
has been used in Ref. || to evaluate the velocity pro- 
file at the surface region in order to evaluate the energy 
loss due to dissipation. In another approach, Ruggeri- 
Friedman [Tl| ] linearize Eq. @ only in the exterior region 
with the assumption p = c 2 e where c s is the velocity of 
sound in the medium. We will discuss about the validity 
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of the above assumption particularly near the station- 
ary configuration as we proceed further. However, the 
advantage of using Eq. (ffl) along with the above relation 
is that one of the variable (say p) can be eliminated so 
that the linear equation becomes a simple wave equa- 
tion. Therefore, we adopt both the approaches as in the 
following. We solve the relativistic hydrodynamic equa- 
tions by linearising Eqs. (|§|) and ([)]) in the interior and 
exterior regions whereas we use the linear Eq. (JlJ) in the 
interfacial region. 



and 

v ( r ) = E. e ~9(r-R) for r > R + ^ (25) 
r 

If k is small, the solution will be the one in which v(r) 
varies slowly over a distance of the order of correlation 
length £ so that g£ << 1. Since k is related to q, next we 
proceed to estimate it by solving the linear hydrodynamic 
equation in the interfacial region and matching it at the 
boundary. 



A. The interior and exterior region 

In these regions, e(r) is varying so slowly that the gra- 
dient energy can be ignored so that Eq. (|ll]) is consistent 
with 



-Vp 



(19) 



Since df/de = at the saddle point [see Eq. (pL3|)] , 
the above relation would imply p is independent of the 
energy density at the stationary configuration. Like 
energy density we also consider a small fluctuation 
in pressure so that we have p(r) — p + (3(r,t) (recall, 
e(r) = e(r) + u(r,t)). We assume that the correspond- 
ing fluctuations satisfy the relation (i = c 2 v where c 2 
is a constant (c s could be the velocity of sound in the 
medium around the saddle configuration). Therefore, in 
the interior and exterior regions, we solve the equation 



K 2 is(r) = c 2 s V 2 is{r) + - I -j] + C ) V a i/(r), 



(20) 



obtained from Eqs. (|^) and (^) after linearizing around 
the stationary configuration and also using the relation 
V 2 /3 = c 2 V 2 ;a Such a relation has also been used in ref 
jnj ) but we differ in our interpretation of the pressure 
gradient. 

Assuming spherically symmetric solutions of the form 



v(r) -- 

we get the relation 
where 



Constant 



-2 2 2 

K = c s q , 



±qr 



R 2 = 



1 



(21) 



(22) 



(23) 



The above relation holds both for QGP and the hadron 
regions, except for the fact that the viscosity coefficients 
are different in two phases. The interior and exterior 
solutions, therefore, are 

v(r) = — sinh(qr) for < r < R - £ (24) 
r 



B. Interfacial region 



As will be shown subsequently, the velocity varies as 
V oc r~ 2 in this region so that r 2 v remains constant. As a 
consequence, V.v = r~ 2 d{r 2 v) / dr vanishes at the surface 
region. Therefore, ignoring the viscous term and elimi- 
nating v from Eqs. (17), an equation for v{r) is obtained 
as 



K 2 v{r) = -V. [Ve (-A-V 2 + /") v(r)\ . 



(26) 



Further, k is assumed to be small so that in the first ap- 
proximation, we can completely neglect the terms con- 
taining K 2 . Thus, to a good approximation in the inter- 
facial region, v(r) satisfies 



dr 



f") u(r) 



= 0. 



(27) 



We follow the procedure described in Ref. to get the 
solution of the above equation in the interfacial region 
(see appendix A for detail), 



v{r) w - 



a(R)R 2 Ae de 



(28) 



where Ae 



2<T r dr 
e^- This solution is quite similar to 
that found in Ref. p] , with n replaced by e and a which 
is now a function of r evaluated at R. Note that in the 
quark and the hadron regions (far away from the surface 
region), e(r) is nearly constant. Thus, Eq. ( p7j ) becomes 
undefined in theses regions. Therefore, a different set of 
equations has been used in the exterior-interior regions 
as discussed in the previous sections. 

We can also get an expression for velocity v(r) from 
the relation, 

kv{t) — — — [r 2 cuv(r)], 



dr 



(29) 



Substituting v(r) from Eq. 

v{r) = — 
r z w 



8J), we get 
de 

rdr—, 
dr 



(30) 



where D is a constant. The above equation can be inte- 
grated to give 



, , ~ D R 

uj q r A 



(31) 



Recall that this result is consistent with our assumption 
that r 2 v is constant in the surface region. 
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C. Dynamical Prefactor 

For the interfacial region, the solution is given in the 
appendix A. It remains now only to apply Eq. ( |l8| ) to 
compute q (or n). As in Q, we can neglect the con- 
tribution coming from the interior region (r < R) and 
the terms of order qR w y^jR in the exterior region. 
The contribution coming from the interfacial region is 
« -aR 3 (Ae) 2 /2a where a is a function of r related to 
the constant B and the second derivative of / w.r.t. e„. 



The exterior region contribution is 
both the terms, we get 



2<jB 



a(R)R 3 (Ae) 2 



B / q . Combining 



(32) 



Assuming a(R + £) = a(R — £) sw a(R) and using the 
relation (<9 2 / /de 2 )^ 1 for B/a [see Eq. (A7)], we obtain 



2a 



q = 



R 3 {Ae) 2 (<9 2 //<9e 2 )' 



(33) 



We can also eliminate d 2 f/de 2 by using the relation 



K de 2 



1 

72' 
"»9 



(34) 



where £ 9 is the correlation length and K is related to the 
surface tension a given by Eq. ( |l2[ ) . The choice of K de- 
pends on the energy density profile e(r) . Following |9] , tr 
can be related to K in the planar interface approximation 
at T r as 



a = K(Ae) 2 /6^ q , 



which will result in 



Q = 



A 

3R 3 ' 



(35) 



(36) 



Therefore, in the case of a non-viscous plasma, we get a 
very simple relation for k given by 



A 

3R 3 



(37) 



where x = £ q /R. This can be viewed as the critical be- 
havior of k that scales as ^~ 1 /(x). However, this scaling 
law is different from the dynamical scaling law that one 
finds in the case of a non-relativistic liquid-vapor tran- 
sition |?],H where k scales as £°R~ 3 . While this needs 
further investigation, one of the reason for this discrep- 
ancy could be unlike the static scaling, the dynamical 
scaling depends on the dynamical behavior of the system 
p6[ which is definitely different depending on whether 
the medium is relativistic or non-relativistic. The above 
result is also valid in the case of the viscous plasma, 
only instead of will scale as £ ? 1 f(x). Therefore, 



for viscous quark gluon plasma, this scaling results in a 
quadratic equation in k with solution given by 



aq 



c s q 



a 2 q 2 
Td 2 ' 



(38) 



where a = (477/3 + £)/u). Since in the first approximation 
g a ft, we can neglect the second term under the square 
root which are higher order in q 2 and a 2 (viscosity). Fi- 
nally, we get 



c s q. 



(39) 



Using Eq. ( pq ) for q, we can obtain a general expression 
for k for a viscous QGP as 



A 

3R 3 



AL 

6R 3 Cu n 



(40) 



Therefore, the prefactor k can be written as the sum of 
two terms having a non- viscous (kq) and a viscous (k v ) 
component. However, both kq and k v have simple de- 
pendence on the correlation length and the bubble 
radius R. As can be seen from Eq. (El), the first term 
is more dominating as compared to the second one par- 
ticularly when T is close to T c . However, as temperature 
decreases, the viscous contribution competes with that 
of the non- viscous one. 

We can also express the above equation (Q) in a dif- 
ferent way by assuming c 2 as 



d 2 f 

Jq d- e r 



(41) 



The above relation is analogous to the non-relativistic 
expression for velocity of sound in the medium which has 
a similar relation with u> and e replaced by n (density) 
O. Then from Eq. (§|) we get 



'2(7 



R 3 (Alo) 2 ' 



(42) 



where we have used the approximation Aw ~ Ae since 
the pressure difference is negligible as compared to the 
difference in energy density. Now using the above q in 
Eq. (^), the prefactor k can be written as 



'2a 



1 



R 3 (Aw) 2 c 2 R 3 (Auj) 2 V3 



(43) 



As can be seen, the first term in the above equation is 
same as Eq. (Q) as obtained by Ruggeri and Friedman 
corresponding to the case of a non-viscous plasma. The 
second term is similar to the result obtained by Csernai 
and Kapusta except with a minor difference, i.e., instead 
of 4, we have a factor of c~ 2 in the numerator [see Eq. (||)] . 
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However this is a small difference which can be removed 
by redefining K [see Eq. (^)J. 

It may be mentioned here that the relation (3 — c 2 s v 
has been used to obtain k as given by Eq. (|4^). This is 
the main result of our work. It is also satisfying to note 
that we can recover the result of Csernai-Kapusta and 
Ruggeri- Friedman under the assumption for c 2 s given by 
Eq. ([hi]). In analogy with the non-relativistic case, we 
may interpret c s as the velocity of sound around the sad- 
dle configuration ( Recall that saddle point is the config- 
uration where Eq. ( p"3| ) is satisfied). Although, the above 
interpretation needs further justification, it is sufficient 
to say that the results of Eq. (^) can be recovered using 
c 2 s as given by Eq. (^l|) . 

D. Result and discussion 

We should point out here that there are several reason- 
able assumptions that have been made in our derivation 
of the dynamical prefactor. As discussed in the text, 
most of these approximations are same as that of used 
in the original work of Langer-Turski || since we use 
the same procedure except that the equations follow rel- 
ativistic hydrodynamics. An important aspect where we 
differ from both Ref || and Ref |6) is the use of the re- 
lation j3 = CgV assumed to be valid in the quark and 
hadron regions, (/3 and v are the radial deviations of the 
pressure and energy density from the stationary solution 
). Although, we differ in our interpretation, such a re- 
lation has also been used by Ruggeri and Friedman ]Tl| ] 
to eliminate one of the hydrodynamic variables. How- 
ever, we do not make any such assumption in the inter- 
facial region. As a result, the linearized equation used 
in the interfacial region is different from the one used in 
the exterior-interior regions. Within above formalism, we 
have derived an expression for the dynamical prefactor k. 
Two important aspects of our result are (a) the prefactor 
k can be written as a linear sum of a non- viscous (kq) and 
a viscous (k v ) component and (b) the non-viscous com- 
ponent (kq) which depends on two parameters R and £ 
is finite in the limit of zero viscosity. The present re- 
sult on kq is also in agreement with the view point of 
Ruggeri-Friedman that the viscosity is not essential for 
the dynamical growth of the hadron bubbles. This fact is 
also evident from Eq. (|^) which implies a non-vanishing 
energy flow wv even in the absence of viscosity. Only 
terms second order in v appear in the energy equation 
in the presence of viscosity and the momentum equation 
contains a term linear in v. This means that viscosity 
terms are relatively unimportant in the energy transport 
for small value of v. The momentum equation, however, 
indicates that viscosity influences the time evolution of 
v. Thus, viscosity can serve to disrupt the energy flow 
and generate entropy but cannot be the only mechanism 
for energy removal. 

The above aspect apparently is in contradiction to the 



general expectation that transport coefficients like vis- 
cosity and thermal conductivity are essential for the re- 
moval of the latent heat and hence for the growth of the 
hadron bubbles []7)-|jo|. This point needs further clari- 
fication. The formation of the hadron bubbles can be 
interpreted as the thermal fluctuation of the new phase 
within a correlated volume of radius R and surface thick- 
ness £. According to the fluctuation-dissipation theory, 
this fluctuation would mean certain amount of heat dis- 
sipation. We can understand the origin of this heat dissi- 
pation as follows. In the absence of thermal conductivity, 
the dissipative losses that occur in a fluid are due to the 
coefficients of shear and bulk viscosity that depend on the 
gradient and divergence of the velocity field respectively. 
In case of an incompressible fluid, the losses are due to 
the shear stress alone since the bulk viscosity that pro- 
vides resistance to the expansion (or contraction) does 
not exist. However, due to the nucleation of the critical 
size hadron bubble, the pressure or the tension in the 
fluid is no longer uniform, the pressure inside the hadron 
bubble being more than the outside. Due to this pres- 
sure difference, the hadron bubbles will keep expanding 
with a non-zero wall velocity. Thus, the fluid medium 
outside will exert a frictional force on the bubble wall 
(causing heat dissipation) whose magnitude depends on 
the pressure difference between the two phases |Q , E8[ . 
In the field theoretical language, this dissipation corre- 
sponds to the coupling of the order parameter <fi to the 
fluid which acts as heat bath. Estimates for it in the 
context of electroweak theory have been given in Refs. 
pj|pO| . Therefore, our non- viscous part of the prefactor 
corresponds to a dissipation of dynamical nature which 
does not depend on any transport coefficients like vis- 
cosity or thermal conductivity. This dissipation basically 
arises due to non-uniform pressure across the interface. 
Following a different approach, Ignatius |H) had also de- 
rived k in the limit of zero fluid velocity to be w 2/(r}Rl) 
where 77 is a phenomenological friction parameter (not to 
be confused with the shear viscosity of the plasma) re- 
sponsible for the energy transportation between the or- 
der parameter and the fluid. Recently, Alamoudi et al 
PH have also studied the dynamical viscosity and the 
growth rate of the nucleating bubble where the viscosity 
effects arise due to the interaction of the unstable co- 
ordinate with the stable fluctuations. They estimate a 
growth rate which depends on R, £ and the self coupling 
A(kk|[1 - .003AT£(f ) 2 ]). In the limit of weak cou- 
pling, the above growth rate scales as which is also 
consistent with our result (23| . 

Finally, we conclude this section with the comment 
that in case of relativistic heavy ion collisions, appre- 
ciable amount of nucleation begins from a super cooled 
metastable QGP phase at which the radius of the criti- 
cal hadron bubble is of the same order as the width of 
the bubble interface. At such point, the homogeneous 
nucleation theory may break down. However, the sys- 
tem moves out of this problematic region quickly due to 



G 



the release of latent heat which heats up the medium 
again towards T c . In the other application, in case of 
cosmology, this problem is not serious, although there 
the actual value of the dynamical growth rate may be 
of less important. In cither case, this study has signifi- 
cance for homogeneous nucleation under the thin-walled 
bubble approximation. 



IV. NUMERICAL RESULTS 

In the following, we compare k obtained from different 
methods. In the case of a second order phase transition, 
the correlation length £ scales in the proximity of the 
critical point as f (T) = f (0)(1 - T /T)- v where v = 0.63 
p4j . However, in the case of a first order phase transi- 
tion, the transition temperature Xb is smaller than T c and 
approaches T c only in the limit when strength of the tran- 
sition becomes weak. Therefore, unlike the second order 
case, £q at T c will be finite and which, in the present 
context, represents the thickness of the interfacial region 
such that R » £ q . Further, we ignore the tempera- 
ture dependence of a and £ g and treat them as constant 
parameters. This assumption can be justified when the 
amount of supercooling is small and the medium returns 
to T c due to the release of latent heat j| . 

Figure 1 shows the temperature dependence of k given 
by Eq. ( ^0|) along with viscous (k v ) and non-viscous (kq) 
components at two different values of a. Following ||, 
we take r\ q as 2.5 T 3 and set ( q to zero. With decreasing 
temperature as well as with decreasing a, the value of 
the critical radius [which is obtained from Laplace for- 
mulae, see Eq. (fhi|)] decreases. Therefore, the k, kq and 
k v increase with decreasing temperature and also they 
have higher values for smaller <r, as expected. The be- 
havior of k v is quite different from that of kq. Initially, 
near T rs T c , the k v has small value, but it exceeds kq as 
temperature comes down particularly at smaller a val- 
ues. Figure 2 shows the similar plot as that of Figure 1 
where we have used Eq. ( |43| ) to estimate k. As seen 
from the figures, both the estimates have similar behav- 
ior although Eq. (|4^)_yields slightly higher values for k 
as compared to Eq.fJhJ). The above studies also suggest 
that the effect of viscosity is negligible at higher a val- 
ues and also for small amount of supercooling. However, 
its effect can not be ignored at much lower temperature 
particularly when a is small. In figure 3. we have also 
compared only the non- viscous part (acq) of the prefac- 
tor as obtained from Eqs.(^) and J43| ) at two different 
values of £ g . Within the present set of parameters, the 
non-viscous parts of the prefactor obtained by both the 
methods behave similar way. 
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FIG. 1. The behavior of K> JUS Qu function of T/T c as ob- 
tained from Eq. (^). Short-dashed curve is for non- viscous 
component kq, long-dashed curve is viscous component k v 
and solid curve is for the sum ko + k v with a — 30 MeV/fm 2 . 
The corresponding curves at a = 10 MeV/fm 2 are shown by 
triangles, squares and circles respectively. 
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FIG. 2. The behavior of function of T/T c as obtained 

from Eq. (M). Lines and symbols are same as in Fig. 1 
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FIG. 3. The behavior of ko (non-viscous component) as a 
function of T/T c . The solid curves are from Eq. ( plo| ) while 
the dashed curve is from Eq. ( [43| ) shown at two typical values 
of£ 9 . 

We also study the dynamics of the nucleation and su- 
per cooling by computing the nucleation rate as 



j K ^0 — Fc/T 



(44) 



where Fq is free energy needed to form a critical bubble 
in the metastable (supercooled) background. The dy- 



namical prefactor k is estimated using Eq. (40) whereas 
the the statistical prefactor Qq is taken from the previous 
works MM as 



fi _ 2 / cr \3/2 / R 
~V~3 V3T/ [j^ 



(45) 



where R is the radius of the critical bubble. Under thin- 
wall approximation Fq and R for a spherical bubble are 
given by 



Fr 



in 



-oR 2 , R 



2(7 



Ph - Pq 



(46) 



From the nucleation rate I(T), the fraction of space 
which has been converted to hadron phase can be cal- 
culated. If the system cools to T c at a proper time r c , 
then at some later time t the fraction h of space which 
has been converted to hadronic gas pj is 



h(r)= / dr'I(T(r'))[l-h(T')}V(T',T). 



(47) 



Here V(t' , r) is the volume of a bubble at time r which 
had been nucleated at an earlier time r'; this takes into 
account the bubble growth. The factor [I — /i(r')] is the 
available space for new bubbles to nucleate. The model 
for bubble growth is simply taken as (25| 



V(T',r) = y(W)) + £dr"v 



(T(r")) 



(48) 



where v(T) = 3[1 -T/T c ] 3 / 2 is the velocity of the bubble 
growth at temperature T Q . (Recall that this velocity is 
different from the velocity of the nucleated bubble surface 
as used in the previous sections). The evolution of the 
energy momentum in 1+1 dimension is given by 



de 



T 



v + C 



(49) 



In this work, we use the bag equation of state for QGP. 
The energy and enthalpy densities in pure QGP and 
hadron phases are taken as 



e q (T) 



3a 9 T 4 



B, uj q (T) 



ia q T' 1 



e h {T) = 3a h T 4 , u h {T) = 4a h T 4 . 



(50) 



(51) 



Here, a q and ah are related to the degrees of freedom op- 
erating in two phases and B is the bag pressure. The 
quark phase is assumed to consist of massless gas of 
?i, d quarks and gluon while the hadron phase contains 
massless pions. Thus the coefficients a q — 377r 2 /90 and 
ah = 37r 2 /90. In the transition region, the energy density 
at a time r can be written in terms of hadronic fraction 
h(r) as 



e(r) = e q (T) + (e h (T) - e q (T))h(r). 



(52) 



Following |15jj3||, the viscosity coefficients for the QGP 
and the hadron phases are chosen as ij q = 2.5T 3 , C, q = 0, 
i]h = 1.5T 3 and (h = T 3 . The other parameters are T c = 
160 MeV and a = 30 MeV/fm 3 . With the above set of 
parameters, Eqs. ([l9]) and (B^) are solved to get h and T 
as a function of time t J^,^ - 

Figures 4(a) and 4(b) show the plot of st- the rate of 
entropy production and T/T c - the rate of supercooling 
as a function of r both for ideal hydrodynamic (IHD) and 
viscous hydrodynamic (VHD) expansions of the system. 
The system cools below T c until nucleation rate becomes 
significant. Afterwards, bubble nucleation and growth 
reheats the system due to the release of latent heat. This 
behavior is similar to what has been studied earlier in 
Refs. [^],[| using a prefactor which explicitly depends on 
the viscosity coefficient of the plasma. In the present 
work, since k has both viscous and non-viscous com- 
ponents, we study the supercooling and extra entropy 
production both with kq and k particularly when the 
medium is non-viscous. First we consider only the ideal 
hydrodynamic expansion. The short-dashed curve and 
the long-dashed curves are obtained using Eq. (^) for 
kq (no viscosity) and k (k v included) respectively. Since 
there is supercooling with kq, extra entropy is generated 
even without viscosity. As shown earlier (see Figures 1 
and 2) , the effect of viscosity on k is not significant with 
a reasonable choice of rj q = 2.5T 3 particularly for small 



8 



amount of supercooling. Therefore, inclusion of k v does 
not effect the supercooling much (see the long-dashed 
curve). The supercooling (hence the entropy production) 
comes down only by about ~ 1% due to viscosity, with 
the present set of parameters. As mentioned before, even- 
though we use K, we do not include viscosity in the hy- 
drodynamical evolution just to bring out the additional 
effect due to the use of k instead of ko in the prefactor. 
However, when the plasma is viscous, the VHD should be 
used for consistency (i.e. when tiy IS included). The use 
of VHD reduces the supercooling by about 10 % as shown 
by the solid curve. Although, the amount of supercool- 
ing reduces, the entropy production goes up. Since the 
effect of viscosity on k is insignificant, the reduction in 
supercooling is purely due to the viscous heating of the 
medium. As a result extra entropy is generated in addi- 
tion to the entropy that is produced due to supercooling. 



b 
b 



(b) 




cal. with fc with IHD 

— cal. with k with IHD- 
cal. with k with VHD 



(a) 



t (fm/c) 

FIG. 4. (a) st- the rate of entropy production and (b) T/T c 
- the rate of supercooling as a function of r, short-dashed 
curves are calculated with fto from Eq. using Ideal Hydro- 
dynamics (IHD), the long-dashed curves are calculated with 
ft using IHD, the solid curves are calculated with ft but using 
viscous Hydrodynamics (VHD). 



V. CONCLUSION 

To summarize, we have derived an expression for the 
dynamical prefactor which governs the initial growth of 
critical size bubbles nucleated in first order phase transi- 
tion. In the case of a non-viscous plasma, the dynamical 
growth rate is found to depend only on the correlation 
length and the size of the hadron bubble which are two 
meaningful scale parameters to describe the critical phe- 
nomena at the transition point. The correction to the 
dynamical prefactor due to viscosity is found to be addi- 
tive and does not affect the growth process significantly 



though additional entropy is generated due to viscous 
heating of the medium. Since the prefactor does not van- 
ish in the limit of zero viscosity, extra entropy is produced 
during the process of nucleation even when the fluid is 
non- viscous. Nearly similar conclusions are also drawn 
by Ruggeri and Friedman who had derived dynamical 
prefactor by solving relativistic hydrodynamics following 
a different approach. However, unlike their result, the 
present prefactor can be written as the sum of viscous 
and non-viscous terms. Interestingly, using an assump- 
tion for velocity of sound in the medium (around the sad- 
dle configuartion) which has a form analogous to what 
is used for non-relativistic plasma, the viscous and the 
non-viscous parts are found to be similar to the results 
as obtained by Csernai-Kapusta and Ruggeri-Friedman 
respectively. 

In the present work we solve relativistic hydrodynamic 
equations both in the interior-exterior, i.e., quark-hadron 
regions and surface regions. The linear hydrodynamic 
equation used in the quark-hadron region is obtained 
after eliminating one of the variable using the relation 
(3 = cr a v which is not valid in the surface region. There- 
fore, a different equation is used for the surface region 
which involves the extra gradient energy. This is where 
we differ from the Csernai and Kapusta method. Further, 
Csernai and Kapusta derived k by equating the flow of 
the outward energy flux with the dissipative loss due to 
viscosity of the medium and the contribution due to a dy- 
namical dissipation was not included. On the other hand, 
Ruggeri and Friedman solve the hydrodynamic equation 
only in the quark region and use a set of boundary con- 
ditions with certain assumptions. In this context, the 
present formalism is more general as we solve the lin- 
earized hydrodynamic equations in all space and obtain 
an expression for the prefactor by matching the solutions 
at the boundary of the interface. Moreover, our result 
is different in the sense that it has a very simple depen- 
dence on the correlation length and radius of the hadron 
bubble although the CK and RF results can be obtained 
from it under certain assumption. 



APPENDIX: A 

We find out the solution for v(r) that satisfies the ra- 
dial equation (see Eq. [27]), 



d_ 

dr 



de 
dr 



+ f" ) X(r) 



= 0, 



(Al) 



where v(r) — x{ r )/ r - Finally we solve for x( r ) from the 
equation 



d 2 d 2 f s 
-*^+fl?U(r)=«(r). 



(A2) 



The above equation is quite identical to the one used 
by Langer-Turski Q in the surface region. The only 
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difference is that the constant a now depends on r as 
a(r) oc (rVe)" 1 . (Note that Ve peaks at r « i?). 
Since i? >> £ and Ve(r) varies sharply in the range 
i? — ^<r<i? + ^, a(r) mostly depends on the Ve(r) 
variation. Therefore, to a good approximation we can 
write a(r) oc (i?Ve) _1 . The general form of the solution 



of Eq. (A2) is given by 



X(r) = / dr'G(r,r')a{r'), 



where G is the Green's function satisfying 



-K- 



dr 2 



de 2 



G{r,r r ) =5{r-r'). 



(A3) 



(A4) 



On either side of the interface d 2 f /de 2 is nearly constant. 
Using the relation V 2 e = 0, it is easy to verify that 



X(r) o(r) 



de 2 



(A5) 



is an approximate solution of Eq. (A2) at the interface 
boundary. Matchi ng the solution in the interfacial re- 



gion given by Eq. (A5) with the solution in the interior 
region given by Eq. (EJ) at R — £ and with the solution in 



the exterior region Eq. (£ 
conditions 



at R + £, give the following 



As'mh(qhR) 



and 



(A6) 



(A7) 



Here the condition q£ « 1 has been used. To get a so- 
lution inside the interface, we follow the same procedure 
as that of Ref. §, i.e., we use the spectral decomposition 
of G as 



Xn(r)Xn(r') 



(A8) 



where A„ are the s-wave eigenvalues and Xn are the cor- 
responding eigenfunctions. For value of r near R, the 
sum will be dominated by the first term. This is be- 
cause Ai ~ — vanishes as R becomes large. Since 
Xi{r) ~ (^■) 1 / 2 (de/dr) is sharply peaked at interface, 
using Eqs. (|A3|) and (|A8|) we get 



X(r) 



a(R)R 2 Ae de 



2a 



dr ' 



(A9) 



where Ac 



e/j. we can also estimate the variation 



of a in the range R — £ < r < R + £, 
a(R + Ve(i?) 



a(R) Ve(R + 



(A10) 



Assuming ^ to be the half width of the full maxima, the 
above ratio could be « y/2. 
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